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Abstract 

Precision measurements of branching fractions and CP asymmetries in non-leptonic 
b — > s decays reveal certain "puzzles" when compared with Standard Model expecta- 
tions based on a global fit of the CKM triangle and general theoretical expectations. 
Without reference to a particular model, we investigate to what extent the (small) 
discrepancies observed in B — > J/tpK, B — > 4>K and B — ► Kir may constrain new 
physics in b — ► sqq operators. In particular, we compare on a quantitative level the 
relative impact of different quark flavours q = c, s, u, d. 



1 Introduction 



Exclusive non-leptonic B meson decays remain a challenge to theory. While semi-leptonic 
B decays are well described within the heavy-mass expansion and allow for a rather precise 
determination of the CKM matrix elements |V^>| an d \V U b\, exclusive non-leptonic decays 
still cannot be described at a similar level of precision. The methods that have been 
proposed so far are based on the flavour symmetries of QCD [1-7], the factorization of 
QCD dynamics in hadronic matrix elements [8-12], or combinations thereof [13-15]. The 
level of precision that one expects from these methods is typically of the order of tens of 
percent, and thus - except for a few "gold-plated" observables - it will in general be hard to 
pin down an effect from new physics (NP) in these decays. Still, from the experimental side, 
the B-factories have collected sufficient information on decay widths and CP asymmetries 
to allow for global fits of the Standard Model (SM) parameters, in particular of the elements 
of the Cabibbo-Kobayashi-Maskawa (CKM) matrix [16,17]. 

The agreement between the standard theory and experimental data is overall satisfac- 
tory, however, in some cases small tensions appear. In the present paper we focus on the 
| AS | = \AS\ = 1 decay modes B — > J/^/K, B —> <pK and B — > Kir, which enter some 
of the present-day "puzzles". Taking the experimental results at face value, we pursue the 
hypothesis that these "discrepancies" with the SM calculations are due to non-standard ef- 
fects [18]. We adopt a model-independent parameterization in terms of isospin amplitudes, 
where we allow for additional contributions from generic NP operators. The moduli, as well 
as the strong phases of the additional terms are then fitted to experimental data on decay 
widths and CP asymmetries. The new weak phase will generally remain undetermined due 
to reparameterization invariance, as long as we do not attempt to fix the hadronic SM 
matrix elements. In the case of B — > Kn decays we make use of additional theoretical 
input from the QCD-improved factorization approach (QCDF) [9]. 

Our paper is organized as follows: In the next section we point out the tensions of the 
SM fit with present data, and give arguments for the way we are going to re-fit the exper- 
imental data including generic NP contributions. In the following section we discuss the 
results of the fits for the individual decay modes and present our conclusions in Section @J 

2 Phenomenology 

2.1 Tensions with the Standard Model? 

In the following we give a brief discussion of the present situation for the B physics ob- 
servables that we are going to consider, where the standard model displays some tension 
with the data: 

• The first point concerns the global fit of the CKM unitarity triangle. Here a small 
mismatch appears between the value of the CKM angle (3 obtained from the direct 
measurement of the time-dependent CP asymmetry in B — > J/i/jKs and the indirect 
determination of the same angle from the mass differences in the neutral B-meson 
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0.758+^ ± 0.075 
(59.6+H ± 5-4)° 



Figure 1: Global fit to CKM parameters from Am^, Am s and |14&/Vd,| ex ci.+inci. • Left: 
Confidence levels in the fj — p plane. Right: Fitted values for CKM parameters, where 
the first error is treated as Gaussian, and the second error is treated as flat. 



systems, Am^/Am^, in combination with the measurement of \V u b/V c b\ from semi- 
leptonic decays [17,19-23]. In fact, using the values from [19], we find for the latter 
case sin 2(3 = 0.758lo o2i ±0.075fl at (see Fig. CD), while the direct determination using 
B J/^K S yields sin 2/3 = 0.678 ± 0.025 [24]. However, the significance of this 
effect depends strongly on the estimates of the theoretical uncertainties, e.g. in the 
determination of |V^&| 5 and can certainly not be taken as a clear evidence for a non- 
standard effect. 

• The second puzzle arises from the time-dependent CP asymmetry in modes like 
B — > 4>Ks which in the SM again yields a determination of sin 2(3, although with less 
precision. The value for (3 obtained from fits to several b — > sss penguin mode£] does 
not agree with the value from the CP asymmetry in B — > J/^Ks [24]. While part 
of the discrepancy may be due to not well understood hadronic effects, it is at least 
curious that the bulk of decay modes involving the b — »■ sss penguin systematically 
yields a lower value for sin 2(3 than the one obtained from B — > J/^/K s (see also [26]). 

• Finally, the theoretical predictions for B — > Kn decay widths and CP asymmetries 
are not always in very good agreement with the data. Within the QCD factoriza- 
tion approach the discrepancy with the data can be brought to an "acceptable" level 
(except for, perhaps, the differences of CP asymmetries AA, see the discussion in 

1 Following the arguments given by HFAG [24], we do not consider the sin 2/3 value extracted from 
B° — > foKs for our discussion, due to the highly non-Gaussian error implied by the BaBar measurement 
[25]. 
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Sec. 13.5.11 below) by assuming particular scenarios within the hadronic parameter 
space, including undetermined l/mj corrections [27]. On the other hand, analyses 
based on SU (3) flavour symmetry for hadronic matrix elements typically have found 
tensions of the order of (2 — 3) a [6,11,13,16,28-31], depending on additional assump- 
tions about hadronic matrix elements. It should be noted that the tensions related 
to the branching fractions have decreased since the inclusion of electromagnetic cor- 
rections in the experimental analysis (for a recent update of the discussion see, for 
instance, [32]). 

Let us, for the moment, take these tensions between theoretical expectations and ex- 
perimental data at face value: Assuming that they are not due to enormous deviations 
from the factorization approximation to hadronic matrix elements, we may try to localize 
in which part of the effective weak Hamiltonian we have to look for NP effects. 

A first possibility are non-standard contributions in the charged b —± u current which 
determines \V U b\. However, it is generally believed to be unlikely that these tree-level pro- 
cesses contain sizeable NP effects. Likewise, the theoretical description of QCD dynamics 
in semi-leptonic decays is fairly well under control, and hence we will not consider this 
possibility here. 

A second explanation could be a non-standard contribution in the mixing phase of the 
\AB\ = 2 part of the effective Hamiltonian, which shifts the observed sin 2(3 to smaller 
values. Such a scenario corresponds to a generalization of Wolfenstein's "superweak inter- 
action" [33]. Obviously, it cannot explain the differences in the sin 2(3 measurements from 
b —> ccs and b — > sss modes. 

The third scenario, which is the one we are going to expand on in this work, is the 
possibility to have an additional contribution in the \AB\ = \AS\ = 1 part of the effective 
Hamiltonian. Evidently, the inclusion of such terms can explain the findings in B — ► 
J/^K S and B —> <t>K$, as well as in B — > Kit. When fitted to experimental data, 
the values for the NP contributions, relative to the leading hadronic amplitudes in the 
SM, can be as large as about 30%. If NP is the explanation for the tensions in non- 
leptonic b —> s transitions, structures beyond minimal-flavour violation (MFV [34-36]) are 
favoured, mainly because the deviations in the B — > (pK CP asymmetries point towards 
an independent NP phase, but also because the constraints on contributions from different 
flavours q in b — > sqq generally can be rather different in size. 
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3 Fit of new physics contributions to experimental data 



3.1 \AB\ = \AS\ = 1 transitions 

Using the unitarity of the CKM matrix, the SM operator basis for non-leptonic b 
transitions can be written as [37] 



rr|AB| = |AS| = l 
-"SM 




c w °S + E c * °< ) + % v ^ v us oJS + E c * °) 

i>3 / \ i>3 / 



(1) 

where 0[ q l are the current-current operators, O^-q the strong penguins, Oj-iq the elec- 
troweak penguins, and CXJ, Cf the electromagnetic and chromomagnetic operators, re- 
spectively. At low energies, the effect of NP in \AB\ = \AS\ = 1 transitions will be 
parameterized by new dimension-six operators. In the following we shall focus on generic 
four-quark operators of the type b —> sqq with q = (6), c, s, u, d, where the Dirac and colour 
structure will not be specified. In order to quantify the possible size of NP contributions, 
we will always assume the dominance of one particular flavour structure, and parameterize 
the corresponding correction to the SM decay amplitudes in a model-independent way. 



3.2 Statistical framework 

The parameter space for the NP amplitudes is explored using the CKMfitter package 
[38]. Here the amplitude parameters are treated as fundamental theoretical quantities, 
and the statistical analysis provides the relative likelihood for a given point in parameter 
space (corresponding to model-dependent "metrology" in the CKMfitter jargon). Other 
theoretical parameters, like hadronic uncertainties from SM physics, are treated using the 
Rfit-scheme, where the corresponding ^-contribution is set to zero within a "theoretically 
acceptable" range, and set to infinity outside. We will sometimes apply the same approach 
to implement additional theoretical constraints/assumptions on the amplitude parameters, 
in order to suppress "unphysical" solutions. 



3.3 Analysis of B -> J/^K 

For B — » J/^K decays the contribution of the second term in the weak effective Hamilto- 
nian (pQ) is small because of two effects: 

• Cabibbo suppression: IV^&V^* |/|Vd,V^,* | ~ A 2 <C 1 

• Penguin suppression: (i) The operators O^. \ do not contain charm quarks, and the 

hadronic matrix elements (J/tpK\0^\B) are suppressed, (ii) The coefficients of the 
loop-induced penguin operators Cj>3 are small with respect to the tree coefficients 
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Table 1: Partial widths [21] and CP asymmetries [24] for B -> J/^K. 



Vcp Sj/i>K s 


-0.678 ±0.026 


Cj/^Ks 

A CP (J/i;K-) 


0.012 ±0.020 
0.015 ±0.017 


T(B~ -> J/-9K-) 
T(B° -> J/M°) 


(6.13 ±0.22) • lO^ps -1 
(5.71 ±0.22) • 10- 4 ps _1 



Furthermore, the electroweak penguin operators can be neglected compared to the strong 
penguin operators. Consequently, in the SM the B — > J/tpK decay amplitude is expected 
to be completely dominated by 

r 6 

MB - J/m = -| V cb V; s (J/^K\C lj2 Og + £ cM e) \B) (2) 

V i=3 

where B = {B®, £>~}, and we projected out the leading [sbcc] component in every oper- 
ator, Oi — > Of\ In particular, the amplitude is dominated by a single weak phase, and 
consequently the time-dependent CP asymmetry in B° — > J/ijiKs is completely deter- 
mined by the B° — B° mixing amplitude, involving the CKM angle (3. Corrections from 
the sub-leading operators have been estimated by perturbative methods at the 6-quark 
scaled and found to give effects of the order of 10~ 3 , only [39,40]. Long-distance penguin 
contributions have been estimated on the basis of experimental data to be not larger than 
10~ 2 [41]. Furthermore, the dominating operators in the SM decay amplitude conserve 
strong isospin (A/ = 0), and therefore do not induce differences between the charged and 
neutral B decays into J/ ijjK. The present experimental data is summarized in Table [D 
We note that the central value for Sj/^k s differs from the indirect determination for sin 2(3 
in Fig. [U but the two values are consistent within the errors. The discrepancy becomes 
slightly more pronounced, if one takes into account the inclusive measurement for \V u b\ 
only, which gives 

sin2/3 = 0.82li°;°*j ± 0.068 flat (using |K&|mci. from [19]). 

Similarly, the central values for the observed isospin-breaking in the CP asymmetries 
{Cj/tpKs vs - — Aqp(J /i!jK~)) and partial widths differ from zero. 

Allowing for generic NP contributions with one weak phase 0\y, the amplitudes can be 
written as 

A(B~ -»■ J/ipK-) = A (B -> J/tfjK) [1 + r e i6w e^° - n e i6w e i(t>1 ] , 
A(B d -> J/ipR°) = Ao(B -> Jji)K) [1 + r e iew e^° + r x e i6w e i(t>1 ] , (3) 

2 The authors of [39] only considered the effect of 0\ 2 . In [40] important contributions from C7 3 _ 6 have 
been included as well. 
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where we have separated the contributions to transitions with AI = (i.e. tree-level 
matrix elements with b — > sec operators or long-distance strong penguins with b — > s(uu + 
dd) or b — > sss operators) and AI = 1 (annihilation topologies with b — > suu or b — > 
sdd). We introduced the absolute values r , r\ and strong phases 0o, 0i for the hadronic 
matrix elements associated with the corresponding NP operators, relative to the leading 
SM amplitude. 



3.3.1 Fit with AI = only (new physics in b — > sec) 

Among the AI = and AI = 1 operators we expect the b — > sec term to give the 
dominating contributions to B — > J/ipK decays, because it has (unsuppressed) tree-level 
matrix elements with the hadronic final state. Therefore, let us first assume that b — > sec 
gives the only relevant NP contribution in which amounts to setting n to zero, while Tq 
should be of the order m^/A^p. Then, the isospin breaking between charged and neutral 
B decays is not affected, and should not be part of the fit. We are thus left with the 
time-dependent CP asymmetries, defined as in [16] 

A CP (f,t) := BR[E °^ /] 1 ( /|~ BR 5^^ / i^ ) := -C f cos(Amt) + S f sm(Amt) (4) 
U ' ' BRffi -> f](t) + BR[B° -> f](t) i V ' } V ' W 

and the direct CP asymmetry A^p(B~ — > J/ijjK~) = —Cj/^k s - Including the contribution 
from r in ([3j), we obtain 



C JH , Ks = -Af v {B~^J/ijK-) 
2r sin O sin 6^ 



1 + 2r cos O cos 6"vk + "^o 



(5) 



2r sin 6^ (cos (2/3) cos <p + r cos(2/3 - 6> w )) 

VcpSj/^Ks = -sm(2/3) + — — —J . 6 

1 + 2r cos 0o cos 6 W + 

We expect the NP amplitudes to provide small corrections to the SM, < r < 1, and 
thus to first approximation we have 

Cj/^Ks — 2r sin W sin O , 
?7cp Sj/^Ks + sin(2/3) ~ 2r sin 6^ cos0 o cos(2/3) . (7) 

From this we read off the interesting parameter combinations 



i . a , y/(vcp Sj/4,k s + sin 2/3) 2 + {Cj /4 , Ks cos2/3) 2 

|r sm^|- _ , (8) 



determining the overall size of the deviations of C from 0, and of S from sin 2/3, and 

Cj/ipK s cos 2/3 



(9) 
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Figure 2: Illustration of the reparameterization invariance: The result for ro smOw 
(left) and tan</>o (right) for the fit to J/ipK observables with NP contributions to AI = 
as a function of cos Ow- (The case of a SM-like NP phase is given by the central values 
cos9\v = —0.38, ro sm.6w = 0.053, tan^o = —0.03, corresponding to the fit in the last 
row of Table [2] below.) 



determining the relative size of the two effects. 

Notice that from the CP asymmetries alone we cannot draw any conclusion about the 
value of the NP phase Ow- This is a consequence of a reparameterization invariance (see 
e.g. [42]) which leaves the decay amplitudes for the neutral B decays in (j3J), as well as the 
branching fraction and the direct CP asymmetry for the charged B decays invariant, 



A' = A(l+£(r e^+r ie ^)) , 

, cos W - £ , sin o w 

cos W = . , sin W = - 

y/l - 2f cosO w + £ 2 V 1 ' 2 ^ cos8 w + £ 



(10) 



and similar transformations for the amplitude parameters r ,i and 00,1, where the parameter 
£ (and therefore also the values for r ,i, 0o,i and Ow) is arbitrary as long as the hadronic 
matrix element Aq for the leading SM contribution is not given explicitly. In particular, 
for r*i = 0, r 1 and small reparameterizations ( < 1, we approximately have 

r - r (1 cos 8 w + 0(e)) , <Po - 0o (l + O^ 2 )) , 
sm6 w -> sin 6 W (l + f cos6 w + 0^ 2 )) , (11) 

which explicitly shows the reparameterization invariance of The reparameterization 
invariance is illustrated in Fig. [2], where as an example we consider the fit result for a 
AI = NP contribution with 6 W = tt — 7sm found in the last row of Table [2] below, and 
apply the reparameterizations in (fTOl) to generate the equivalent solutions for other values of 
6w- In particular, we verify that the combinations r sin Ow and tan0 o are approximately 
reparameterization-invariant, except for Ow near zero or 7T. 

As a consequence of the reparameterization invariance, the fit to the experimental 
data will generally allow for "unphysical" solutions, where the strong and weak phases are 
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tuned in such a way that the absolute size of the NP contribution r can be unreasonably 
large. In order to suppress such effects, we implement additional constraints: (i) For small 
NP contributions, the fit should not depend on the parameter combination |r cos 9w\; 
constraining |r cos 6^1 < 0.4 should therefore only affect the unphysical solutions, (ii) If 
the phase 9w of the NP operator is close to the SM one, we do not expect to be sensitive 
to NP in CP asymmetries in any case; we may therefore concentrate on 30° < 9w < 150°. 
(iii) For 9w = ^ — 7sm our fit can also be interpreted as a determination of the size of 
sub-leading SM contributions from Cabibbo- and penguin-suppressed amplitudes, which 
possibly may have been underestimated in [39,40]. In this case, one could also include 
the information from B — > J/tpir decays to further constrain the hadronic parameters, 
using SU(3) flavour symmetry [41], and correcting for the different relative CKM factors. 
Considering the CP asymmetries in B — > J /ipn alone, we find that the constraints on r e l5 ° 
are less restrictive than and consistent with the B —> J/ipK case. The ratio of branching 
fractions in B — > J/ipir and B — > J/i/jK further constrains r [41]. However, we find that 
this ratio essentially depends on the combination 

1 T[B° -> J/jjK } „ A 2 

2 r[B°-J/^r°] ~ R 2 sm \± + rr 1 ' 

and thus the constraints on r are highly correlated with the assumptions on the 577(3) 
breaking parameter Rsu{3) f° r th e ratio of the leading B — > J/i[)ir and B — > JjijjK am- 
plitudes. As this ratio cannot be estimated in a model-independent way at present, we 
refrain from a detailed quantitative analysis. However, it should be mentioned that for 
Rsu(3) ~ l s smaller values for r are favoured. 

Using the experimental values for Cj/^Kg, Sj/^Ks, and Aqp(B~ — > J/ipK~), together 
with the value for sin 2(3 from the indirect determination in Fig. [U we fit the preferred 
ranges for the NP parameters - applying the different constraints as discussed above - 
as shown in Table [2] and Fig. [3] Since the value for \V u b\ from the average of inclusive 
and exclusive decays is close to its indirect determination from sin 2/3, the fitted range for 
r sm9w in this case is consistent with zero, and the related strong phase O is uncon- 
strained. Still, for sufficiently small strong phases, NP contributions of the order 20% are 
not excluded either. On the other hand, taking into account the inclusive value of \V u b\ 
(with its small tension with sin 2/3) only, the fit prefers non-zero values for r sm9w of the 
order 5-30% and relatively small strong phases 0o- (Notice that small strong phases are 
generally expected within the QCD factorization approach to hadronic matrix elements 
in the heavy-quark limit [8].) Compared to the estimate of SM corrections in [39,40], 
the typical order of magnitude for r is thus significantly larger. Although the present 
experimental situation is not conclusive, our analysis shows that an improvement of the 
experimental precision for B — > J/ ipK observables or of the theoretical precision in the 
\V u b/V c b\ determination may still lead to interesting conclusions. 
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Figure 3: Fit results <po vs. ro sinf?^/ for different scenarios, see also Table [2j The plots 
on the upper half refer to the case where \ V u b\ is determined from exclusive and inclusive 
decays, whereas for the plots in the lower half only the inclusive value is used. In the 
plots on the left only the constraint |ro cos 9w\ < 0.4 is imposed. The plots on the right 
are for fixed values 9\y = vr — 7sm ■ 
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Table 2: Fit to direct and mixing-induced CP asymmetries in B — ► J/ipK, using the 
indirect determination of sin 2/3 and including the AI = NP contribution ro, only. We 
show the lcr confidence level for the two relevant parameter combinations \ro sin 6w\ and 
00, using different additional constraints to suppress "unphysical" solutions (see text). 
The upper half of the table corresponds to using the sin 2(3 value from the indirect fit 
with IVu&lexci.+inci. in Fig. [TJ In the lower half, only |Vub|inci. from [19] is used. 



Scenario 


r sin^ 


tan 0o 


excl.+incl. 9w free \tqCOs9w 
30°<9 W <150° \r cos9 w 
30° < B w < 150° \r cos9 w 
W = tt-7sm \r cos9 w 


< 0.4 
free 

< 0.4 
free 


[0 to 0.23] 
[0 to 0.19] 
[0 to 0.19] 
[0 to 0.13] 


unconstrained 
unconstrained 
unconstrained 
unconstrained 


incl. 6w free r^cos^ 
30°<^<150° \r cos9 w 
30°<^<150° \r cos6 w 
6> w = 7r-7sM \r Q cos9 w 


< 0.4 
free 

< 0.4 
free 


[0.02 to 0.34] 
[0.03 to 0.33] 
[0.03 to 0.33] 
[0.03 to 0.19] 


[-0.41 to 0.18] 
[-0.26 to 0.12] 
[-0.26 to 0.12] 
[-0.24 to 0.11] 



3.3.2 Fit with AI = 0, 1 (new physics in b —> suu or b — > sdd) 

New physics contributions to either b — > suu or b — > sdd may lead to isospin asymmetries 
between charged and neutral B — > J/ipK decay rates and CP asymmetries. In this case 
we may fit with both r ^ and T\ ^ 0, and consider the observables in Fig. [3] together 
with the (CP-averaged) isospin breaking in the decay rates [21] 

MB J/rpK) - r[Bd - J j- j ^- ] - r[B± - J j—^ ] - -0.035 ±0.026. (13) 

For small values of r and r±, following [43], we have the approximate relations 

Vcp S + sin 2/3 ~ 2 (r cos0o + T\ cos0i) sin9w cos 2/3 , 

A^p ~ — 2 r sin O sin Q w , 

AAcp — —2 n sin 0i sin 9w , 

Aj ~ 2 ri cos 0i cos 9 W . (14) 

They are manifestly invariant under the approximate reparameterizations, following from 
(TT0|) in the limit f = 0(r o ,i) < 1, 

sinflw -> sin# w (l +£ cos9 w + C>(£ 2 )) , 

cos 9 W — > cos 0iy — £ sm2 #w + C(£ 2 ) , 

r cos 0o + r\ cos 0i — > (r cos0o + r x cos0i) (l - ^ cos^vy + C(£ 2 )) , 

r\ cos 0i — > r\ cos 0i (l + £ sin 6* w tan 9\y + C(£ 2 )) , 

r ,i sin 0o,i -> r ,i sin0 o ,i (l - £ cos0 w + C(£ 2 )) . (15) 
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Figure 4: The result for tq e l< ^° (left) and n e i<?il (right) in the complex plane from the 
fit to J /il)K observables, with isospin-breaking NP contributions 6 — ► su-u or 6 — > sdd. 
The new weak phase has been fixed to cj)w = vr — 7sm- 



To keep the discussion simple, we may again concentrate on the special case 9\y — tt — ^y. 
The fit result is plotted in Fig. The la parameter ranges are given by 

r cos 0o = [-0.077 to 0.112] , r sin0 o = [-0.008 to 0.006] , 

r lC os0! = [ 0.013 to 0.088] , n sin^ = [ 0.000 to 0.015] . (16) 

Notice that again, the strong phases for the preferred ranges turn out to be small. Solutions 
for other values of 6w can be reconstructed by means of the reparameterization invariance 

We conclude that small deviations from the SM expectations in B —> J/ipK can be 
explained by NP in either b — > suu or b —> sdd, alone. However, one has to keep in 
mind that, compared to the contributions from b — > sec, the b — > suu or b — >• sdd only 
contribute via penguin (ro) or annihilation (7*1) diagrams to hadronic matrix elements. 
Thus, an additional suppression with respect to the tree-level matrix elements fitted in the 
last section (see Table [2]) is expected. Notice that, depending on the actual size of these 
suppression factors, our result for r and r\ may also be interpreted as due to unexpectedly 
large effects from sub-leading SM operators. Again, the information from B — > J/iJm 
observables together with assumptions on SU(3) breaking effects could be used to further 
constrain r n and r\ in this case. 

3.4 Analysis of B — ► qbK 

The discussion of B — > (j)K decays is very similar to the B — > J/ipK case. The most 
important difference is due to the fact that a tree-level operator for b — > sss transitions is 
absent in the SM, and therefore the leading SM amplitude Ao(B — > <pK) already receives a 
penguin suppression factor of order A compared to Aq(B — > J/ijjK) (see for instance [44]). 
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Consequently, the relative size of both, Cabibbo suppressed SM contributions as well as 
potential NP contributions, may be enhanced accordingly. Indeed, the experimentally 
observed discrepancy between S^k s an d sin 2/3 is more pronounced, while estimates within 
the SM typically give small effects [5,12,45-47]. 

To keep the notation simple, we use the same symbols r i? 0j as in the B — > J/ipK to 
parameterize NP contributions to the B — > <pK decay amplitudes 

A(B -> <j)K) = A (B -> <f>K) [1 + r e i9w e i(t>0 T n e^e** 1 ] . (17) 

However, one has to keep in mind that both, the involved NP operators and the strong 
dynamics in hadronic matrix elements, are different. 

3.4.1 Fit with AI = (new physics in b — > sss) 

Using the experimental values for the direct and mixing-induced CP asymmetries in B — ► 
4>K together with the value for sin 2(3 from the indirect determination in Fig. [TJ we fit the 
preferred ranges for the NP parameters as shown in Fig. [51 Again, we only quote the result 
for a particular value for the new weak phase, 9w = it — 7sm- Other solutions follow from 
the same reparameterization invariance as in (fTOl) . Comparison with the B — > J/ipK case 
in Fig. [3] shows: 

• Again, the fit prefers small strong phases (j> . 

• The preferred value for tq in B — > <pK is by a factor of 2-3 larger than the one in 
B — > J/ipK. After correcting for the penguin suppression factor, phase space and 
normalization, this implies that the coefficients of the involved NP operators in both 
cases may be of similar size. 

We emphasize, that the latter observation also implies that unusually large hadronic pen- 
guin matrix elements in the SM could simultaneously explain the B — > J /ipK and B — > (pK 
discrepancies. 

3.4.2 Including AI = 1 operators 

The current data shows no evidence for isospin asymmetries in B — > (f>K decays [24], 

AA CP (B^(pK) = 0.02 ±0.13, 

Ai{B^4)K) = 0.04 ±0.08, (18) 

although again the relative effects from b — > suu and b — > sdd operators are expected to 
be larger than in the B — > J/ipK case. We find it instructive to turn the argument around 
and estimate the potential size of isospin violation in B — > (pK by simply rescaling the 
solutions for r and r\ in (JTBi) by a factor 2.5 (see above), which yields the "1-cr estimates" 

AA CP (B -> <f>K) ~ (0 to 0.14) , 
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quantity 


value 


0.6 


sin 2(3 


0.758t£^ ± 0.075 


0.5 


VCP S<pK s 


-0.39 ±0.18 


0.4 




0.01 ±0.13 


0.3 




0.034 ± 0.044 




\r singly 


[ 0.08 to 0.35 ] (la) 


0.2 


tan 0o 


[-0.24 to 0.11 ] (Iff) 



ro sin Qw 

Figure 5: Fit to direct and mixing-induced CP asymmetries in B — > <f>K, using the 
indirect determination of sin 2(3 and including the contribution of a NP operator with 
AI = 0, i.e. b — ► sss. The NP weak phase is set to 9\y = vr — 7sm- Left: Confidence 
levels for the two relevant parameter combinations | rosing I and 0o- Right: Input 
parameters (upper half [24]) and Iff ranges for the output values (lower half) of the fit. 



A^B — ► (bK) ~ -(0.17 to 0.01) . (19) 

The resulting order of magnitude is comparable with the present experimental uncertain- 
ties. If our estimate makes sense, a moderate improvement of the experimental sensitivity 
could already lead to a positive signal for isospin violation in B — * (bK. 

3.5 Analysis of B — » Ktt 

In the SM, the general isospin decomposition for B —> Kn decays can be parameterized 
as [4,9] 

A(B- -> 7T~K°) = P(l + e a e i(f>a e"* 7 ) , 
-V2A(B~ -> ir°K-) = P(l + e a e 1 ^ e"* 7 - e 3/2 (e~ 17 - qe iuJ )) , 

-A{B d -> n + K ~) = P(l + e a e^e- i7 -e r e^ (e~^ - q c e iuJc )) (20) 

and 

V2A(B d -> 7T°/? ) = A(B~ -> n~K°) + V^-A^" -> 7T°/i-) - A(5 d -> TT + if") 

fixed by isospin symmetry (i.e. neglecting QED and light quark-mass corrections in the 
hadronic matrix elements). Here P is the dominating penguin amplitude, whereas the 



13 



quantities €7,3/2 contain tree-operators but are doubly CKM-suppressed. Without any 
assumptions on strong interaction dynamics, in the isospin limit one is left with 11 inde- 
pendent hadronic parameters for 9 observables. In order to test the SM against possible 
NP effects in these decays, one needs additional dynamical input. Qualitative results from 
QCDF [9] include: 



The SU(3)f symmetry prediction [3] 



only receives small corrections. 



3 \v cb v;\ c 9 + c 
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2 \v ub v*\ d + C 2 



(21) 



The parameter e a e z ^ a is negligible in QCDF. Consequently the direct CP asymmetry 
in B~ — > tt~K° is tiny (in accord with experiment). 

The parameter qc e luJC is of minor numerical importance. 

The parameters and e 3 / 2 are expected to be of the order 20-30%, with the related 
strong phases of the order 10°. Furthermore, at least at NLO accuracy, the differ- 
ence between e^ T and €3/2 e 4 ^ 3 / 2 is a sub-leading effect proportional to the small 
coefficients a 2 7,9 in QCDF. 



In the subsequent fits, we will set e a to zero and use the values from [9], 



q = 0.59 ± 0.12 ± 0.07 , lu = -0.044 ± 0.049 , 
q c = 0.083 ± 0.017 ± 0.045 , to c = -1-05 ± 0.86 , 



(22) 
(23) 



in order to reduce the number of independent hadronic parameters within the SM to 5. 
(Notice that the overall penguin amplitude parameter P in (|20|) will not be constrained 



from theory, but will essentially be fixed by the experimental data for the B z 



7T 



branching fractions.) Tensions in the fit, or incompatible values for the parameters eT, 3 / 2 
and 0t,3/2 then may be taken as indication for possible NP contributions. 



3.5.1 New physics in B — > Kw ? 

The critical observables in B — » Ktt transitions are [32] 

"BR(_B~ -> ir°K~) + BR(B + -> vr ^) 
BR(B- -> tt'K ) + BR(B+ -»■ n+K°) 
BR(B d -> ix+K-) + BR{B d -> tx-R + ) 



R c 

Rn 

AA 

C n o Ks 



2 [ BR(B d -> 7T°K ) + BR(B d -> 7r ir ) 



^cp(-^ - - ^cp(^ - Tf^K r 

0.14 ±0.11, VcpS 7T o Ks = -0.38 ±0.19 



1.11 ±0.07. 



0.97 ±0.07. 



0.142 ±0.029 



(24) 
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Table 3: SM fit results for ex, 4>t, e 3 / 2 , 4> 3 / 2 , with e a = and qe luJ and qc e %u)c varied 
according to (I22|23|) from [9]. The best fit values for the latter parameter are obtained as 
q = 0.49, ijj = 0.005, qc = 0.038, uc = —1.91. For comparison, we show in the last line 
estimates for the corresponding hadronic parameters from [32] which have been obtained 
by relating B — > ttK to B — > irir via SU (3) relations and dynamical assumptions (central 
values only). 





e T 


4>T 


£3/2 


03/2 


ReAe 


Im Ae 


Best: 


0.21 


0.21 


0.04 


0.07 


0.18 


0.07 


la: 


[0.10, 0.32] 


[0.10,0.50] 


[0.01,0.15] 


[0.05,0.09] 


[0.07,0.33] 


[0.05,0.09] 


2a: 


[0.05,0.44] 


[0.05,1.32] 


[0.00,0.38] 


[0.03,0.11] 


[-0.13,0.42] 


[0.03,0.11] 




r 


5 


r c 




-p n cos6» n 


-p n sin^ 


[32] 


0.12 


0.44 


0.20 


0.02 


-0.10 


0.04 



Within our SM approximation, we expect (see also [13]) 

R c -R n ~ 2e 3/2 (e T -e 3/2 (l-q 2 ))+0(\ 3 ), (25) 
AA ~ C n o Ks ~ 2 (e T sin T - e 3 / 2 sin 3/2 ) + C(A 3 ) , (26) 
77cpSW s - -sin 2/3 + 2 cos 2/5 (e T - e 3/2 ) + 0(A 2 ) , (27) 

where we used that 6^,3/2 ~ A, 0r,3/2 ~ A, g c ~ 0, u; ~ 0, and cos 7 ~ A in the SM. 
Considering the recent experimental data, the first relation turns out to be well fulfilled, 
whereas the second and third relation require a sizeable difference between ere^ and 

To quantify this observation, we perform a fit (within the SM) to the quantities e-r e l ^ T 
and e 3 / 2 e"^ 3 / 2 , as shown in Table [3l In Tabled] (3rd column) we compare the best fit result 
with experimental data and observe a very good agreement. In particular, the expected 
approximate equality AA ~ C 7r o Ks is fulfilled by the data. The fitted values for the 
individual amplitude parameters 6t, 4>t, £3/2, 03/2 are in qualitative agreement with the 
expectations from QCDF. However, the comparison of e T e l< ^ T and e 3 / 2 e^ 3 / 2 shows sizeable 
deviations, 

Ae:=e T e^ T -e 3/2 e^ ^0, 

which are incompatible with the NLO predictions from QCDF (for the status of NNLO 
predictions, see [48-51]). In the notation for topological amplitudes [52] this would corre- 
spond tclH a ratio C/T = |Ae/er| in the range [0.52 — 3.00] with the central value at 0.89. 
Assuming that higher-order QCD effects and non-factorizable power corrections cannot 
substantially change the approximate equality between and e 3 / 2 , this might be taken 
as a weak indication of NP in B — > Kn decays (for a recent discussion, see also [31]). It 

3 In a parameterization where the annihilation topology A is explicit [4], one has |Ae/e<r| = (C+A)/(T— 
A), where T{C) denote the colour-allowed (-suppressed) tree amplitude. 
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Re(Ae) 

Figure 6: SM fit results for Ae in the complex plane, with e a = and qe lLU and qc e lu ° 
taken from [9]. 



is also interesting to compare the fitted values for Ae with the latest estimates obtained 
in [32] on the basis of SU(3) relations and dynamical assumptions about sub-leading decay 
topologies, see last row in Table [3] In this case, a sizeable C/T ratio is obtained from a 
fit to the B — > nix observables, but with the "wrong" sign for the corresponding strong 
amplitude, compared to our SM fit. As the dynamical mechanism for generating (sizeable) 
strong phases in charmless non-leptonic B decays is not completely understood, a resolu- 
tion of the observed discrepancies in Ae from non-factorizable QCD corrections within the 
SM cannot be excluded (see, for instance, the discussion in [53]). 

We may interpret the required difference between e^ and e 3 / 2 as due to NP contributions 
in the AJ = 1 Hamiltonian. In this case the fit result for the quantity Ae, shown in Fig. [6], 
is already a measure for the possible effect of NP operators. Notice however, that again 
the weak phase associated with these operators cannot be fixed. To continue, we follow a 
similar line as in the analysis of B — > J/ijjK and B — > <f)K decays, and assume that only 
one particular NP operator of the type b — > sqq gives a significant contribution in B —> Kn 
decays. 

3.5.2 New physics contributions with AI = only 

The presence of a NP contribution with AI = (in our case, this includes the "charm 
penguin" b —> sec, as well as b — > sss and b — > s{uu + dd)) has the same impact as the SM 
parameter e a in (|20|) . except for a possibly different weak phase. Within our approximation, 
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Table 4: Experimental data for B — > i^7r-decays vs. various best fit results. The third 
column shows the SM fit with Ae / 0, which corresponds to x 2 /d.o.f. = 2.43/3. The 
fourth column shows the best fit result for Ae = (with ex e** T = e3/2 et ^ 3/2 varied 
according to their QCDF ranges, see text) and a NP contribution with AJ = and 
6w = vr — 7sm, yielding x 2 /d.o.f. = 18.5/6. The last column shows the analogous fit 
result for a NP contribution from (essentially) b — ► suu with 6w = vr — 7sm, which 
corresponds to x 2 /d.o.f. = 2.91/3. Experimental values taken from HFAG [24]. 



Observable 


HFAG (after ICHEP'06) 


SM fit 


NP (7 = 0) 


NP (7 = 0,1) 


BR(tt u K-) • 10 b 


12.8 ±0.6 


12.2 


12.6 


12.6 


BR(tt-K ) ■ 10 6 


23.1 ± 1.0 


23.9 


23.8 


23.8 


BR(ir + K-) ■ 10 6 


19.4 ±0.6 


19.7 


19.6 


19.6 


BR(7r°K°) ■ 10 6 


10.0 ±0.6 


9.5 


9.0 


9.2 


Acp(7t-K») 


0.009 ± 0.025 


0* 


-0.02 


0* 


A C p(7r°K-) 


0.047 ±0.026 


0.048 


0.001 


0.049 


Acp{-k+K-) 


-0.095 ±0.015 


-0.095 


-0.06 


-0.094 




-0.38 ±0.19 


-0.39 


-0.34 


-0.48 


C n o Ks 


0.12 ±0.11 


0.14 


0.06 


0.13 


R c 


1.11 ±0.07 


1.02 


1.06 


1.06 


Rn 


0.97 ±0.07 


1.04 


1.09 


1.07 


AA 


0.142 ±0.029 


0.143 


0.06 


0.143 



one thus obtains 

?(l+r e*e % ) , 

P (1 + r e Ww - e 3/2 (e~^ - qe w )) , 
P (1 ± roe** e idw - e T e^ T (e" i7 - q c e iulc )) , (28) 

where r e l ^°e l9w parameterizes the NP amplitude with A7 = 0. As explained above, the 
QCDF approach predicts small values Ae ~ 0. In the following NP fits to B — > Kn decays, 
we will therefore fix Ae = for simplicity, and vary the common values in the ranges 

e T = e 3/2 = 0.23 ± 0.06 flat ± 0.05 gauss , <f) T = <f) 3/2 = -0.13 ± 0.1 l flat , (29) 

which have been determined by combining the QCDF errors [9] on the individual parame- 
ters (flat errors are combined linearly, and the larger of the Gaussian errors is chosen). As in 
the B —> <j)K example, since the leading SM amplitudes are already penguin-suppressed, 
we expect r < 0(1) and (f> < 0(A). Generically, we now expect a sizeable direct CP 
asymmetry in B~ — > tt~K° of the order A. The experimental value for that asymmetry 
should therefore be included in the fit and will essentially constrain the parameter combi- 
nation r sin 0o- On the other hand, using the power-counting ei,qc,u,4>i ~ A, a A7 = 
NP operator does not contribute to the critical observables Aj and AAcp in (1241) at order 



A{B~ -> 7T"7?°) ~ 
-V2A(B~ -> 7T°7\") ~ 
-A(B d -> 7T+A'") ~ 
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Table 5: Fit to AI = NP contribution in B — > Kir. We show the la confidence levels, 
assuming Ae = (with er e l ^ T = e 3 / 2 e 1 ^/ 2 varied according to their QCDF ranges, see 
text), e a = and with qe tuJ and qc e luJc varied according to (|22|23p from [9]. 



9w 


l r o| 


tan 0o 


X 2 /d.o.f. 


57r/6 


[0.31 to 0.43] 


[0.00 to 0.03] 


14.9/6 


2vr/3 


[0.23 to 0.35] 


[0.01 to 0.06] 


17.9/6 


7T - 7SM 


[0.22 to 0.34] 


[0.01 to 0.07] 


18.5/6 


tt/3 


[0.23 to 0.50] 


[0.06 to 0.15] 


24.6/6 


7r/6 


[0.15 to 0.68] 


[0.21 to 0.54] 


34.4/6 



A, either. As explained in [32] and references therein, these observables are sensitive to 
AI = 1 operators which, in the SM, are represented by electroweak penguin operators. 

As a result, the NP fit with AI = contributions generally leads to a bad description 
of the experimental data, except for certain fine-tuned parameter combinational with small 
NP phase 6 W and unreasonably large values for the amplitude normalization factor P. To 
avoid such fine-tuned scenarios, we consider some particular examples with fixed NP phase 
9 W , see Table [5l We thus confirm on a quantitative level that AI = NP contributions 
alone cannot resolve the B —> Kn "puzzles". 



3.5.3 New physics with AI = 0, 1 (b — > suu or b — > sdd) 

New physics contributions with AI = 1 induce two new isospin amplitudes 

P, 



r <W e *w e i4F*> P f and r (3/2) e *w e ^ 3/2; 



corresponding to final \Ktt) state with I — 1/2 or 1 = 3/2. Using the connection between 
( 1201) and isospin amplitudes (see e.g. [53]), we obtain (again within our approximation) 



A(B~ 
■V2A(B- 



7l~K°) 
7T°K-) 



P I 



J<t>0 



(3/2) 



D i9w 



P(l+r e i9w - e 3/2 (e~^ - q e iu} ) 



-A{B d -> ix + K~) ~ P l + r e 



100 a i6w 



sm e < /2) + r (3/2) ei0 (3 m 



J8v 



(30) 



In order to reduce the number of free parameters in the fit, and to avoid unphysical 
solutions, we apply additional assumptions/approximations: 



4 Notice, that contrary to the B — > <j>K and B — ► J/ij)K analyses, we cannot exploit reparameteri- 
zation invariance here, because we decided to constrain certain hadronic input values from QCDF. As a 
consequence, the fit results will explicitly depend on the value of the NP weak phase. 
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Table 6: Same as Table H for the fit with AI = 0, 1 NP contribution in B -> i^7r. 





| r (l/2)| 


tan0i 1/2) 


| r (3/2), 
I'l 


tan<^ /2) 


xVd-o.f. 


5tt/6 


[0.04 to 0.08] 


[0.06 to 0.08] 


[0.00 to 0.04] 


unconstr. 


4.3/3 


2tt/3 


[0.03 to 0.07] 


[-2.65 to -0.51] 


[0.00 to 0.05] 


unconstr. 


3.5/3 


7T - 7SM 


[0.03 to 0.09] 


[-9.89 to 0.38] 


[0.00 to 0.07] 


unconstr. 


2.9/3 


tt/3 


[0.04 to 0.11] 


[-16.4 to -0.42] 


[0.41 to 0.51] 


unconstr. 


0.4/3 


7f/6 


[0.20 to 0.26] 


[-0.41 to -0.03] 


[0.65 to 0.70] 


unconstr. 


1.7/3 



• Following the experimental observation, we force the direct CP asymmetry in B — > 
7T~K° to vanish identically, which yields the relation 

r e^+rf /2) + rf /2) = , 

which we use to eliminate the parameters r and 0o- This effectively implies that we 
deal with a b — > suu operator which does not contribute to B~ — > n~K° in the naive 
factorization approximation. 

• Again, we assume the SM contributions to the amplitude parameters and e 3 / 2 to 
lie within the QCDF ranges, see (1291 ). 

(1/2) -A 1 / 2 ) (3/2) va' 3 / 2 ) 

In Fig. [7J we display the results for the NP parameters r[ e t(pl and r\ e t(pl in 
the complex plane, for different values of the NP weak phase $w- The corresponding la 
ranges are collected in Table O The resulting central values for the observables in the case 
0w = 7T — 7sm are listed in the last column of Table EH We observe that the fit depends 
on the value of the NP weak phase 9w in an essential way. In particular, depending on 
whether 6w is less or greater than 7r/2, we encounter disjunct regions in parameter space. 
One of the regions always corresponds to relatively small values of r [ 1 ^ 2,3 ^ 2 ' ) < 10%, whereas 
for values of 9w close to or tt solutions with r ^ 2 ' 3 ^ as large 50% are possible. 
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1 /2 'j, I S/2 'Apt 

Figure 7: Fit results for AI = 1 NP contributions r/ e^i (left) and jy e 4<?, i (right), 
with e a = 0, Ae = and Aqp(it~ K°) = 0, see also text. The plots in the upper row 
refer to a weak phase 6\y = 5ir/6, the ones in the middle row to 6\y = n — 7, and the 
lower ones to 0w = tt/6. 
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4 Conclusions 



To date, flavour physics is evolving from the -B-factory era to the LHC era. While the 
former has led to an enormously successful confirmation of the CKM mechanism in the 
SM, the latter is expected to reveal direct and indirect signs for physics beyond the SM 
with interesting interplay between high-p^ and flavour physics [55-57]. In this context, a 
crucial task is to constrain the flavour structure of NP models, manifesting itself in rare 
quark and lepton decays and production and decay of new flavoured particles. 

While within concrete NP models the chiral, flavour and colour structure of new op- 
erators could be completely specified, the present work pursues a model-independent ap- 
proach. Assuming the dominance of an individual NP operator, the analysis of B — > J/t^K, 
B — > (f>K and B — > Kn observables allows us to infer semi-quantitative information about 
the relative size of NP contributions to b — > s cc, b — »• s ss, b — ► s dd, and b —> suu 
operators. The main conclusions to be drawn are: 

• From the comparison of isospin-averaged B — > J/ipK and B — > (pK observables 
we find that - after correcting for relative penguin, phase-space and normalization 
factors - NP contributions to b — > sec and b — > sss operators may be of similar size 
(order 10% relative to a SM tree operator). 

• In a scenario, where b — > sdd is the only source for NP contributions in B — ► ttK 
observables, while the SM contributions are estimated in QCD factorization, one 
cannot simultaneously explain the individual CP asymmetries. In particular, the 
experimental value for Acp(vr + K~), which does not receive leading NP contributions 
from b — > sdd, cannot be reproduced in a scenario with negative strong phase 4>t- 

Moreover, the small direct CP asymmetry for B~ — > tt~K° requires the matrix 
element of a b — > sdd NP operator to have either a small coefficient or a small phase. 

• This leaves the b — > suu operators, which correlate isospin-violating observables in 
B —> J/ipK and B — > Kn decays, and may be even somewhat larger (order 20% 
relative to a SM tree operator) than the b — > sec and b — > sss operators. 

In all cases, in order to explain deviations from SM expectations for CP asymmetries 
without fine-tuning of hadronic parameters (see the discussion after (TTTT) ) . we have to 
require non-trivial weak phases (0w 7^ 0, 7r), which could be due to NP, albeit the case 
0w = 7T — 7sm is always allowed, too. Consequently, our findings are still compatible with 
a SM scenario where non-factorizable QCD dynamics in matrix elements of sub-leading 
operators is unexpectedly large. 

In the future, an improvement of experimental accuracy, in particular on the isospin- 
violating observables, could lead to even more interesting constraints on the relative impor- 
tance of different b — > sqq operators and their interpretation within particular NP models 
with MFV [34-36] or beyond (see e.g. [58-62]). 
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